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ABSTRACT 

We analytically derive the correlation functions of the stochastic response 
of period-p orbits in a generic one-dimensional map system under the influ- 
ence of weak external, additive white noise. Two approaches, a recurrence 
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scheme and a path integral method, are provided. Beside the supercycle, we 
recognise that the fluctuations on each elements of the period-p orbits are 



colored noise which have non-vanishing time correlation. The results also 
indicate that the stochastic responses will be stationary after a large num- 
ber of period-p iterations. The analytical results are confirmed by numerical 
simulations. 
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1 Introduction 



The role of dynamical noise in the deterministic nonlinear system has been 
studied extensively [Q - • These results indicate that, in a noisy dynamical 
system, not only the higher periods become obliterated by the increasing 
noise level but the bifurcation points themselves become blurred. In partic- 
ular, at the points of bifurcation, the Lyapunov exponents no longer vanish 
as they do in the deterministic limit. Meanwhile, through both numerical 
and analytical studies, some universal scaling behaviors are discovered at 
the onset of chaos even in the noisy dynamical system @-||. Mayer- Kress 
and Haken carried out numerical investigation on the distribution function 
of stochastic response of a logistic model with uniform white noise using an 
integral recursion formula || . But the analytical form of the correlation func- 
tion, or distribution function, of the stochastic response of period-p orbits in 
a generic noisy map system was not presented. 

Later on, in a review paper jTj, Crutchfield, Farmer and Huberman 
studied the stochastic response of a noisy one-dimensional map, x n+ i = 
rf(x n ) + p n , where f(x) = x(l — x) and p n is the external noise. They 
proposed that the stochastic response should be regarded as the fluctuation 
on the control parameter r. It means that the noisy map can be rewritten as 
x n+ i = (r + q n )f(x n ), where q n is the response parametric noise. However, 
when we numerically study the bifurcation diagrams of the noisy logistic map 
and the noisy Ikeda map, see Fig. 1, we notice that the mean values of x n on 
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each elements of noisy stable orbits coincide with those of the stable orbits 
without noise. This fact implies that the dispersion should be regarded as the 
fluctuation on the elements of noise-free period-p orbits {x$ | 1 < m < p} 
of the one-dimensional map x n+ i = f r (x n ). In this map, the element of 
period-p orbits x$ and control parameter r obey the periodicity condition, 
x$ = and the stability condition, \M&\ = \ Ul=i(^\^))\ < 1. 

Based on this observation, we will derive the general formula for the 
correlation function of each elements on the period-p orbits by a recurrence 
scheme and by a path integral method. The results indicate that, beside 
the supercycles, the stochastic responses are colored in nature which means 
that the correlation between the kp-th step and the (k + £)p-th step are non- 
vanishing. It also shows that the correlation functions, or the distribution 
functions, of the stochastic response are stationary after kp iterations when 
k » 1. The plan of this paper is as following. We derive the correlation 
functions by recurrence scheme in section 2, and by path integral method in 
section 3. Numerical verification on two interesting systems, logistic map and 
Ikeda map, are demonstrated in section 4. Finally, some concluding remarks 
are provided. 
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2 Recurrence scheme 



Here, we consider a generic one-dimensional nonlinear map in which the 
external additive noise is involved dynamically, 

Xn+l — fr(Xn) + Vn, (1) 

where r\ n is white noise with zero mean, i.e. < i] n T] m >= DS nm and < rj n >= 

2 

0. For examples, D is ^- for uniform white noise with noise level a and 
D = a 2 for Gaussian white noise with standard deviation a. 

To derive the stochastic response on each elements of fluctuating period-p 
orbit, we let the initial state x m to be located in the neighbourhood of x$ , 

1. e. x m = x$ + £ m , £ m << I. After kp iterations, the stochastic response on 
the m-th element of period-p can be regarded as the fluctuation on the m-th 
element of noise-free period-p orbits, that is 

Xkp+m = %kp+m £,kp+m = + ^kp+m- (2) 

Here, ^kp+m is the stochastic response on the m-th elements of period-p after 
k times of period-p iteration. 

From now on, the noise level and stochastic response will be assumed 
to be small in comparison with the noise-free period-p orbits such that the 
perturbation scheme is applicable. By substituting Eq.(2) and X(k+i)p+ m = 
%rn + £(fc+i>+m into the p times iteration of Eq.(l), i.e. 

%(k+l)p+m fr\---fr\fr\%kp+m) Vkp+m ) + Vkp+m+1---) + Vkp+m+p-1- (3) 
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It leads to a recurrence relation between ^ k+ i) p+m and £ kp + m , 

t {k+l)p+m = M^i kp+m + H% m 1) . (4) 
Here, we have denoted 

a^^SjA (5) 



i=Q 



and 



Vk+m , forp = l, 

E^=i 1 (Mj^ 1) ^ fcp+m+i _ 1 ) + 7/ fej)+m+ j ) _i , for p > 2, 



where 



and 



ji" s n*. <8) 



By iterating the recurrence relation, Eq.(4), one obtains 



a kP+ m = {M^tu + Y.(M^r i Ht^. ( 9 ) 

i=l 



For simplicity, but without loss of generality, we set the initial state to be 



•Em "^m 



that is £ m = 0. Therefore, the stochastic response, after kp times 



of iteration, becomes 

i kp+m = E(mW)*- 1 ^- ( 10 ) 
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As a result, the mean fluctuation of stochastic response is zero when the 
< Vk >= 0. Here, < rjk >= limAr_ >00 -^ J^fjo 1 Vk ^ denotes the ensemble 
average. The correlation function of stochastic response between the kp-ih 
step and the {k + £)p-th step can be written as 

i'=l i=l 

(11) 

Using the properties of the input white noise, we find, for the case p > 2, 

^ 11 {k+t)-i' ,m 11 k-i,m 
p— 1 p— 1 

= E ^2 M f,m' ^jfm^ < V(k+e~i')p+m+j'-lV(k-i)p+m+j-l > 
j'=lj=l 

+ < r l{k+Z-i')p+m+p-l r l(k-i)p+m+p-l 

= E E M&^^SJ'kv-tD + ^hi'-e-D 
j'=ij=i 

= (l+g(M& 1) ) 2 )^'^- ( 12 ) 

Therefore, the correlation function ( for p > 2 ) between the fcp-th step and 
the (A; + £)p-th step is 

<f p , (^ (p) ) £ (i - (^ (p) ) 2fc ) fl | ^V^-D^n n ^ 

It should be noted that the term E P jZl(M^~ 1] ) 2 in Eq.(13) will be zero for 
the most trivial case p — 1. 
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3 Path integral approach 

As it has been shown in previous section using the recurrence scheme, 
one can derive the correlation function of stochastic response of a dynamical 
map contaminated by uniform or Gaussian white noise. Due to the excellent 
form of Gaussian white noise 

1 rj 2 

piVi) = v^ exp ^^ (14) 

one also can derive its stochastic response by using path integral approach. 

For simplicity, we assume the initial state is located at x m = x^\ it 
means that the probability distribution of x is P(x m ) = 5(x m — x$) at the 
beginning. After one iteration, x m+ i = f r (x m )+Vm, the distribution function 
of x m+ \ disturbed by the Gaussian noise has become 

P(Xm+l) = J dx m p(x m+1 - f r (x m ))P(x m ) (15) 
1 , ( Xm+1 - f r (x$))^ 

= v^ exp ( ^ }■ 

After the next iteration, x m+2 = f r ( x m+i) + Vm+i, on e can write down the 
distribution function of x m+ 2 

P(x m+2 ) = J dx m J dx m+1 p(x m+1 -f r \x m ))p(x m+2 -fr(x m+1 ))P(x m ). (16) 

Now, we introduce the kernel K(x m+ e, x m ) which is defined as 

P{x m+e ) = J dx m K(x m+e ,x m )P(x m ). (17) 



The kernel can be read out from Eq.(16) and (17), that is 

-l e-i 



K(x m+ e,x m ) — / dx m+ j Y[ K(x m+ j + i, x m+ j), (18) 



j=l j=0 



where 



K(x m+J+1 ,x m+3 ) = -^exp{- ^± J^^ 2 }. (ig) 
Then, the distribution function of Xk p + m , after kp iterations begins with 



(20) 

To obtain the distribution function by integration, we change the variable 
x m+ j to x m+ j = x^ + j + £ m+J '. As mentioned in section 2, £ m+ j is the cor- 
responding stochastic response and is assumed to be small enough that the 
linear perturbation scheme is applicable. Therefore, we have dx m+ j = d£ m+ j 
and 



fr{x m +j) — fr( x m+j) + o ( P ) £m+j + 0(£ ) 



dx 



m+j 



x 



m+j+1 



+ (21) 



and the distribution function can be expressed as 

(22) 
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Using the formula 

1 r (x-az) 2 } 1 f (z-byf 



2a. 



2 



1 exvi ( x ~ ab y) 2 ) (23) 



1 _ (M^) 2fc - - 



and integrating out d£ m ...d£(kp-i)+m step by step. Finally, we have 

P{x kp+ m) = /o ^ 2 exp{- {Xkp+ 2 y^ ))2 }, (24) 

where 

<»» = A^Jiw){ 1 + £<*£">'}■ < 25 > 

This result indicates that the stochastic response on each elements of period- 
p orbit for the input Gaussian noise still is a Gaussian-like distribution with 
mean x$ and deviation o"fe p+m . This is consistent with the case i = of 
Eq.(13). 

It is expected that the stochastic response will inherit the time correlation 
from the deterministic dynamical system. Here, we define the time correla- 
tion between the kp-th step and the (k + £)p-th step of stochastic responses 
as 

< { x (k+i)p+m~ < x (k+t)p+m >){ x kp+m~ < x kp+m >) > 
= < £,(k+£)p+m£,kp+m > 

J d,X (k+£)p+m J ^^fcp+m£(fc+^)p+m£fcp+m^(-£(A;+^)p+mj x kp+m) P {, x kp+m) ■ 

(26) 



From Eq.(18) and Eq.(23), one finds 

K(x {k+ e) p+m , x kp+m ) = - exp\ * — '- \, (27) 

where 

Therefore, the correlation between the (A; + £)p-th step and the kp-th step is 

< i { kU) P+ mtkp +m >= (M^Ya 2 kp+m . (29) 

These results are exactly the same as Eq.(13) which is obtained from the 
recurrence scheme. 

4 Numerical verification 

To check the correctness of our analytical derivation, we compare the 
analytical results to the numerical results in two interesting systems: logistic 
map and Ikeda map. Firstly, let us consider the noisy logistic map which 
is contaminated by uniform white noise with amplitude a = 0.001. Fig. 2(a) 
shows that the ratio between the standard deviation of stochastic response 
and input noise, y / <l ^ fcp+TT ^ fcp+m> , on each m-th element of period-1, 2, 4 orbits 
in the noisy logistic map, x n +i = rx n {l — x n ) + i] n , for 1.0 < r < 3.54. The 
normalized correlation function between the (k + £)p-th step and the kp-th 
step ^ ^ <g (fc +0P+^ P +^> ; on each m _ th e i ement Q f per iod-2 orbit in the logistic 

map for £ = 1,2, are presented in Fig. 2(b) and 2(c), respectively. Similar 

10 



comparison for the case of the noisy Ikeda map, x n+ ± = 1 — r(l — cos(jX n )) + 
r] n , are demonstrated in Fig.3(a)-(c). 

When the Gaussian white noise is added into the dynamical system, it 
is worth examining the distribution function P(xk p + m ) directly. Fig. 4 (a)- (c) 
show the numerical and analytical distribution functions P(xk P + m ) of ra- 
th element of period-2 orbit in the logistic model with control parameters 
r = 3.1, 1 + \/5 and 3.3, respectively. 

As one can see from the graphs, all the numerical results agree with the 
analytical derivation very well. So long as the ratio between the amplitude of 
the input noise and the value of fixed point is less than 5%, and the dynamical 
system is not too close to the bifurcation point. 

5 Concluding remarks 

From the derivation of the general formula for the correlation function 
of each elements on the period-p orbits by the recurrence scheme and the 
path integral method, we know that the recurrence approach is good for 
any white noise with zero mean, but path integral approach is only feasible 
for Gaussian white noise. However, the path integral method give us the 
distribution function, with this, one can easily to write down the higher 
moments of correlation function directly. Beside the supercycles (M^ = 0), 
the stochastic response on each elements of period-p orbit are colored which 
have the correlation time rjf 1 = ~~ ln i^(p)r This time correlation can be 
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regarded as an inheritance from the deterministic properties of the dynamical 
system. 

The general formula Eq.(13) also indicates that the amplitude of stochas- 
tic response is always larger than the level of input noise . The response on 
the supercycle is the most attractive one that the response of each element 
is independent of k, the times of period-p iteration, that is 

< >= {1 + J2( M j p m 1] ) 2 } D - (so) 

3=1 

The minimum response occurs at the m = 1 elements on the supercycles in 
which M^ l) = and the response is identical to the original input white 
noise, see Fig. 2(a) and Fig.3(a). 

Moreover, from the stability condition \M^\ < 1, one can deduce that 

< Z kp+ U( k+ e) P+ m >^ 1 _ {M { P)) 2 { 1 + TXM^y}D, (31) 

for k » 1. It means that after a large number of period-p iterations the 
correlation functions on each elements become stationary. Therefore, when 
one defines the long-time average after large k (transient) as 

^ N-l 

< £,kp+m£,kp+m >LT= j^^Tt E £(k+i)p+m£,(k+i)p+m ■ (32) 

i=0 

The correlation function almost matches that computed by the ensemble 
average. Fig.5 shows that Eq.(30) is a good approximation for the correlation 
function evaluated by the long-time average of the stochastic response of 
logistic map after the transient steps k = 1,000 , if the neighborhood of 
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the bifurcation point is excluded. This approximation has been reported in 



our previous papers [|T(J JTT|. Indeed, those numerical results also indicate 
that the noise leads to a shift of the bifurcation points. This feature could 
account for the failure of our schemes near the bifurcation points. It may 
imply that one should regard the response of noise as parametric noise in 
the neighborhoods of the bifurcation points, such as Crutchfield had done 
before, or deal with the response by other method |T^| . 

As an application, we have used the recurrence scheme to study the 
stochastic response of colered noise and the noise cascading in uni-directionally 



coupled elements system ||11|| . Those results offer us an opportunity to de- 
velop suitable methods to control noise in the system, see reference [11]. 

The recurrence method can be extended to study the stochastic response 
of periodic orbit in continuous time system, as long as one can write down the 
Poincare map of the dynamical system. In fact, Svensmark and Samuelsen JT2] 
had obtained a similar result for period- 1 case in a noisy dynamical system, 
see Eq.(ll) in reference [12]. They also found other type of auto-correlation 
function in a noisy dynamical system under the influence of a period-2 res- 
onant perturbation. On the other hand, using the Floquet theory, Wiesen- 
feld JL3] obtained the noise response of a special model of continuous time 
system as an integral of noise source and fundamental matrix, and give a 
generic correlation function for the period- 1 orbit. Both papers focus on the 
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influence of noise in the vicinity of period-doubling, but pay little attention 
on the correlation function for each elements of period-p orbits. 

Note in addition: After the completion of this manuscript, Prof. Weiss 
informed us his related works published about a decade ago |14| . In his paper, 
he calculated the first and the second moments of the probability distribution 
for a N-dimensional noisy map, up to the second order approximation. The 
i = case in our results is a special case, the first order approximation of 
one-dimensional noisy map, of his result in slightly different form, but he 
did not state the response of each elements on the stable orbits explicitly. 
Moreover, he did not compute the time correlation function between the kp- 
th step and (k + £)p-th step for the case £ ^ 0, that reveals the inheritance 
of the dynamical property from the deterministic structure. 
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Figure Captions 



Figure 1. 

The coincidence between the mean value of each elements of noisy stable 
orbit and those stable orbits without noise, are demonstrated in (a) noisy 
logistic map, x n+ i = rx n (l — x n ) + r/ n , and (b) noisy Ikeda map, x n+ i = 
1 — r(l — cos(^x n )) + r/ n . The amplitude of white noise r\ n is 0.001. The cross 
(x) denotes the average for 10,000 ensembles on each elements of period- 
p orbit after 1,000 times of period-p iteration. The solid line denotes the 
elements of noise-free orbit . 
Figure 2. 



(c), respectively, see text. Here, the dynamical system is the noisy logistic 
map with noise amplitude a = 0.001. The diamond (o) denotes the numerical 
simulation for 10,000 ensemble and the solid lines are the analytical results. 
Figure 3. 

Similar to Fig. 2 (a)-(c), but the Ikeda map x n+ i = 1— r(l — cos(^x n ))+r) n , 
go instead of the logistic map. 
Figure 4. 

The distribution function P(xk p + m ) of the ra-th element of period-2 in 
logistic map with control parameters (a)r = 3.1, (b)l + a/5 and (c) r = 3.3. 
Here, k = 10 and the standard deviation of input Gaussian white noise is 
a = 0.1. The diamond (o) denotes the numerical simulation and the solid 




via r for k = 10, £ = 



,1,2 are shown in (a), (b) and 
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lines are analytical results. 
Figure 5 

The diamond (o) denotes the numerical results of correlation function 
evaluated by long-time average for different r in logistic map. Here, we take 
50,000 data for the long-time average after 1,000 transient steps. The solid 
lines are the analytical results of Eq.(30) for k — \, 000, £ = 0. 
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